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Abstract: We study the metal abundances of F and G type stars in the Galactic disk and halo using 
the Canada— France— Hawaii Telescope Legacy Survey D4 field. For a sample of stars within 7 kpc 
of the Galactic plane, we derive mean abundance values of [Fe/H]= —0.77 ± 0.36 dex for the thick 
disk, and [Fe/H]= —1.42 ± 0.98 dex for the stellar halo, respectively. These metallicites are consistent 
with the metallicity estimates from the Sloan Digital Sky Survey Data Release 6. We do not find a 
vertical abundance gradient for the thick disk between 1 kpc and 4 kpc. However, stars within 1 kpc of 
the Galactic plane are saturated in our dataset, and therefore we cannot rule out the observed vertical 
metallicity gradient for the thick disk stars in the SDSS. Moreover, we find a negative trend in the halo 
metallicity with increasing distance from the Galactic plane. This trend could be due to a contribution 
from an increasing number of very metal poor stars that belong to the outer halo. However, systematic 
effects in photometric metallicities are largest for the most metal poor stars, and therefore these effects 
cannot be ruled out. 
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1 Introduction 

Our view of the formation of the Milky Way has changed 
dramatically since the discoveries of complex substruc- 
tures in the Galaxy. The Milky Way has four spatially 
and kinematically distinct components; thin disk, thick 
disk, halo, and the bulge. However, recent discoveries 
of irregular structures like the Sagitarrius dwarf tidal 
stream and the Monoceros stream show that the for- 
mation of galaxies is not a steady process resulting 
in a smooth distribution of stars. Instead, galaxies are 
constantly shaped by the infalling smaller galaxies (see 
Juric et al. 2008 and the references there in). 

The Sloan Digital Sky Survey (SDSS) and the 
Sloan Extension for Galactic Understanding and Ex- 
ploration (SEGUE) are providing photometry and spec- 
troscopy for a very large number of stars that can be 
used to study the stellar populations in the Galaxy. 
Such a study involving ~ 20, 000 stars with the SDSS 
spectra demonstrated that the Galactic halo is likely 
to have two components; an inner and an outer com- 
ponent (Carollo et al. 2007). These two components 
seem to have different spa tial density profiles, s tellar 
orbits, and metallicities. ICarollo et al. I (|2007h sug- 
gest that the outer halo formed through dissipation- 
less chaotic merging of smaller subsystems within a 
pre-existing dark matter halo , supporting the com plex 
galaxy formation scenarios. Ive zic et al. I (|2008h use 
the SDSS photometry data to study the metallicity 
distribution of disk and halo stars. However, their 
dataset is not deep enough to test the dichotomy of 
the halo. 

The Canada— France— Hawaii Telescope Legacy Sur- 
vey (CFHTLS) provides a valuable new source to 



study the metallicity distribution of the Galactic stel- 
lar populations. The third data release of the Deep 
survey now provides ugriz photometry down to a lim- 
iting magnitude of 26. CFHTLS lack spectroscopic 
observations, which would provide the most accurate 
metallicities. However, imaging data can also be used 
to deduce stellar metallicities through the traditional 
UV excess method (Wallerstein 1962; Sandage 1969). 
The metal absorption lines mostly affect the UV part 
of the spectrum, therefore metal— rich and metal— poor 
stars with the same effective temperature (or the same 
B — V color) can be differentiated based on their UV 
colors. In addition, imaging data enables us to ob- 
tain a complete flux-limited sample of stars and push 
the studies of stellar metallicities to fainter magni- 
tudes than the spectroscopic surveys. Here we use 
the CFHTLS data for the D4 Field to perform such 
an analysis. Section 2 describes the data and the 
selection of point sources. The derivation of metal- 
licites from the photometric data and its application to 
the CFHTLS ugriz system is discussed in Section 3, 
whereas the results from this analysis are discussed in 
Section 4. 



2 The Data 

The CFHTLS D4 field is located at RA = 22 h 15 m 31 s , 
DEC = -17°43'56" (J2000), I = 39°, and b = -53°. 
The TERAPIX T003 public data release for the D4 
field covers an area of 0.86 deg 2 and provides photom- 
etry with limiting AB magnitudes of u = 26.5, g = 
26.3, r = 26.4, i — 26.0, and z = 25.0, respectively. 
The high galactic latitude of this field is similar to the 
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Figure 1: HLR versus z-band for the D4 field. 
Note that objects with i < 17.00 are saturated. 
Star— galaxy separation seems to be robust down 
to i = 22 mag. The solid horizontal line marks the 
boundary for selection of point sources (HLR < 
2.82). 



fields targeted by the SDSS. 

Without follow-up spectroscopy, we have to rely 
on morphological classificatio n to select point sou rces. 
ISchultheis et al. I (|2006f > and lLimboz et al. I (120081 ) showed 
that the star— galaxy separation is possible through 
the use of the half— light radius (hereafter, HLR) mea- 
surements from the i-band images. HLR corresponds 
to the radius that encloses 50% of the total flux from 
the source. Figure 1 shows the HLR versus i-band 
magnitudes for the sources in the D4 field. Objects 
with i < 17 mag are saturated, and therefore have 
unreliable HLR measurements. Using spectroscopic 
data from the VLT-V IRMOS Deep Survey (VVDS; 
Le Fevre et al. 2005), iLimboz et al~l (120081 ) demon- 
strated that the spectroscopically confirmed stars have 
HLR(i) < 2.82 and that the star— galaxy separation is 
reliable down to i — 21 mag. Since the CFHTLS goes 
five magnitudes deeper than this, we now look into 
pushing the point source classification limit to fainter 
magnitudes. 

Figure 2 presents (u — g) versus (g — r) color— color 
diagrams for four different magnitude bins. The top 
left panel shows the color— color diagram for point sources 
(black points) and resolved sources (red points) for 
i — 20 — 21 mag, where the point source classifica- 
tion seems reliable. This panel reveals a tight stellar 
sequence, and confirms the expectations that stars can 
be identified reliably down to 21 mag. The other pan- 
els in this figure show the color— color diagrams for 
i = 21 - 22, 22 - 23 and 23 - 24 mag. These three 
panels show that the stellar sequence is still clearly 
visible for i = 21 — 22, but the galaxy contamination 
starts to become a problem for i > 22 mag. We count 
the number of sources that are classified as point-like 
and that lie to the left of the stellar sequence, and 
compare it with the number distribution of resolved 
sources with HLR(i) > 2.82 to estimate contamination 



of the stellar sequence by galaxies. We estimate that 
the galaxy contamination is less than 4% in the range 
(g — r) — 0.1 — 0.7 mag for i — 21 — 22 mag. An analysis 
of the number distribution of resolved and un-resolved 
sources in the other panels show that the galaxy con- 
tamination ranges from 7% to 23% for i = 22 — 23 
mag, and even higher for i > 23 mag. Therefore, we 
restric our sample to the sources with i < 22 mag and 
HLR < 2.82. Using the above criteria, we classify 

7348 sources as stars. 

We use the dust maps of ISchlegel at al. I (|1998| ) 
to de-redden the photometry. E(B — V) for the D4 
field ranges from 0.023 to 0.031 mag, with a mean of 
0.027 mag. The small range of reddening correction 
for the D4 field shows that this correction does not 
have a signifincant effect on our analysis. Perhaps the 
most important correction required for the photome- 
try is the correction from the MegaCam system to 
the standard ugriz filter system. Studying the differ- 
ences b etween the u g riz sy stem and the MegaCam 
system, IClem et al~l (|2008l ) find that the MegaCam 
griz filters are consistent with the standard system, 
though the u band data need to be corrected using a 
third order polynomial. The MegaCam u— band fil- 
ter was selected to take advantage of the better UV 
transparency of the CFHT and it is slightly redde r 
than the SDSS u- band filter. IClem et al. 1 (|2008f ) 
compare theoretical isochrones in the standard sys- 
tem and the corrected MegaCam system and they 
find that the isochrones overlap each other within 0.01 
mag for both dwarfs and giants with colors u — g = 
0.0 to +3.0 mag and g — i = —0.5 to +3.5 mag, and 



= -2.3 to 0.0. We use the 
Clem et al. I (|2008D (see their 



with metallicities [Fe/H 
transformation given by 
Figure 5) to convert the MegaCam u— band photom- 
etry into the standard u— band system. In addition, 
a correction is required to transform the photometry 
from the primed system to the SDSS 2.5 m telescope's 
unprime d ugriz system . We use the transformations 
given by iTucker et al. I ||2006D to transform the pho- 
tometry to the unprimed system. 
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Figure 2: (u—g) versus (g — r) color— color diagram 
for different i-filtcr magnitude limits. Black points 
shows stars (HLR < 2.82), whereas red points rep- 
resent galaxies (HLR > 2.82). 
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The apparent g magnitude distribution for our sam- 
ple of stars is shown in Figure 3a. The number of 
stars increases with increasing magnitude until g — 
22.1 mag, which we use as the completeness limit. For 
a ~ G5 main— sequence star with (g—r) = 0.4 mag (see 
Table 1 in Ivezic et al. 2008), this limit corresponds 
to a Galactocentric distance of 21 kpc and a Galactic 
plane distance of «17 kpc. Figure 3b shows the (g — r) 
color distribution for our sample. There are two dis- 
tinct peaks located at ~ 0.20 mag and 1.15 mag, corre- 
sponding to F stars and late type M stars, respectively. 
We select our sample of F and G type stars by restrict- 
ing the sample to the stars with 0.20 < (g — r) < 0.60 
(shown by vertical dashed lines). This results in a 
sample of 1315 F and G stars. Majority of these stars 
lie abov e the s ynthetic (u — g) vs. (g — r) relation of 
iPickles I |l998), indicating that they are metal— poor 
disk and halo stars. 



3 The Metallicity Distribution 



Using SDSS spectroscopy and photometry. Ilvezic et al.l 
(2008) have developed a metallicity estimation method 
based on the (it — g) and (g — r) colors. Their method 
reproduces the metallicities obtained from low reso- 
lution SDSS spectroscopy with a root-mean-square 
scatter of 0.2 dex. They also derive a photometric par- 
allax relation using ugriz photometry of globular clus- 
ters. The agreement between the CFHT /MegaCaM 
and the SDSS ugriz filter sets, after correcting the 
MegaCam it— band photometry (Clem et al. 2008), 
enables us to use their photometric metallicity and 
parallax methods on our dataset. Photometric metal 
abundances of the F and G dwarfs with 0.2 < (g — 
r) < 0.4, are determ ined via the following equation of 
Ilvezic et al.l (|2008l ). 



[Fe/H] = -4.37 - 8.56(u - g) + 15.5(g - r) 
-39.0(u - g)(g - r) + 23.5(u - gf + 20.5(g - r) 2 



+12.1(« - gf(g -r) + 7.33(u - g)(g - 
-W.l(u-g) 3 - 21.4( ff 



(1) 



For G dwarfs with (g — r) > 0.4, we replaced the 
(u — g) color in the above equation with (u — g) — 2(g — 
r) + 0.80. We estimate the absolute magnitudes of 
our sample o f stars using the following relations from 
Ilvezic et al.l (120081 1. 



M r (g - i, [Fe/H]) = M r (g - i) + AM r ([Fe/H]) (2) 



AM r ([Fe/H}) = 4.50 - l.ll[Fe/H] - 0.18[Fe/H} 2 (3) 



M r (g 



-5.06 + 14.32( 



12.97( fl 



+6.127(p - if - 1.267( 5 - if + 0.0967( 5 - i) 5 (4) 

Figure 4 displays the u — g color distribution and 
the derived metallicities for our sample of F and G 
stars. Both of these distributions are bim odal, and 
they are fairly well fitted with two gaussians. Ilvezic et al. I 
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Figure 3: (a) g magnitude distribution. Drop in 
g = 22.10 is taken as the completeness limit, (b) 
(g — r) distribution. The presence of the two main 
peaks at (g — r) ~ 0.20 and ~ 1.15 are seen. 
Vertical dashed lines show the color index range 
0.20 < (g — r) < 0.60, which corresponds to F and 
G type stars. 



(|2008D also found a bimodal u — g color distribution for 
the SDSS sample with two peaks at u — g m 0.9 and 
1.1 mag. The two peaks observed in the color distri- 
bution of our sample are located at u — g — 0.86 ±0.13 
and 1.10 ± 0.41 mag. The agreement between the col- 
ors of our sample of F and G stars and that of the 
SDSS F and G stars shows that the transformations 
between the MegaCam photometry and the standard 
ugriz system are reliable. The two peaks in the u — g 
color distribution correspond to metal— poor halo stars 
and more metal— rich and closer disk stars, respec- 
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Figure 4: (a) (u — g) color and (b) [Fe/H] distri- 
bution for our sample of 1315 F and G type stars. 
Solid lines represent two Gaussian distributions fit- 
ted to the data. 



tively. The width of 0.13 mag is significantly larger 
than the median u — g error of 0.02 mag at g — 22 
mag, providing a measure of the intrinsic variations in 
the u — g color and the metallicity distribution. 

The observed peaks in the u — g color distribution 
translates into two peaks at [Fe/H] = —0.78 ±0.39 and 
— 1.62 ± 0.82 in the metallicity distribution. Figure 5 
shows this distribution as a function of increasing dis- 
tance from the Galactic plane. The mean metal abun- 
dances for the fitted gaussians are also listed in Table 
1. The mean metal abundance shifts as a function of 
galactic plane distance. The thick disk dominates at 
z — 1 — 4 kpc, and a small contribution is still visible 
up to 8 — 12 kpc, after which only the halo stars with 
[Fe/H] w -1.7 are visible. 
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Figure 5: The Metallicity distribution as a func- 
tion of the distance from the plane of the Galaxy. 



4 Discussion 

Our mean metallicity estimate of [Fe/H] = — 0.78±0.39 
for the thick disk is consistent with the results from 
the analysis of F and G stars in the SPSS. Using the 
SPSS Data Release 3 spectroscopy. I Allende Prieto et al~1 
(2006) found that thick disk G dwarfs with 1 < z < 3 
kpc show a compact metallicity distribution with a 
peak at [Fe/H] ~ —0.70. Our mean metallicity esti- 
mates for 2 = 1—2, 2—3, and 3—4 kpc (see Table 1) 
range from —0.71 to —0.73 dex, and they are entirely 
consistent with the Allende Prieto et al. result. How- 
ever, the mean metallicity for the thick disk is more 
metal— poor compared t o the literat ure values of —0.60 
(|Gilmore feWvse Il985h and -0.50 (|Carnev et al. 1 19891 ). 
Since the SPSS photometric metallicity relation is 
known to underestimate the metallicity at high metal- 
licities (see the discussion in Ivezic et al. 2008), our 
mean metallicity estimate for the thick disk probably 
suffers from these systematic effects. 

U sing the SPSS Data Release 6 photometry. Ilvezic et al.l 
(|2008l ) found that the median metallicity for the thick 
disk smoothly decreases with distance from the plane 
from [Fe/H] = -0.60 at 500 pc to -0.80 at several 
kpc. We do not observe such a metallicity gradient 
for the thick disk between z = 1 kpc and 4 kpc, but 
our mean metallicity estimate is consistent with their 
mean metallicity estimate beyond several kpc, —0.80 
dex. In fact, if we restrict our sample to z < 7 kpc, 
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Figure 6: The Mean metallicity for the thick disk 
(filled dots) and halo (open triangles) as a function 
of distance from the Galactic plane. 



then the mean metallicity peak for the thick disk stars 
in our sample shifts to [Fe/H] = —0.77 ±0.36, entirely 
consistent with the Ivezic et al. analysis. Since stars 
with i < 17 mag are saturated in the CFHTLS data, 
we do not have any stars within 1 kpc of the Galac- 
tic plane. The trend found by Ivezic et al. is mostly 
based on the relatively brighter stars in their samplqj, 
which show a peak in metallicity at [Fe/H] ~ —0.60. 
The majority of the z > 1 kpc thick disk stars in their 
sample show a peak around [Fe/H] ~ —0.70 (see their 
Figure 7), similar to our analysis. 

Figure 6 shows the metallicity trends for the thick 
disk and halo as a function of the Galactic plane dis- 
tance. If we ignore the thick disk metallicity estimates 
beyond 8 kpc due to small number statistics, then a 
small shift in metallicity is observed for the thick disk 
between 4 and 8 kpc. However, the mean metallicities 
for the stars with z = 1—4 kpc and 4—8 kpc are consis- 
tent within the errors. In addition, we do not find any 
significant trend between the Galactocentric distance 
and metallicity for the thick disk stars. 

Our metallicity estimate of — 1.62T0.82 for the halo 
i s consistent with the literature value s of —1.60 dex 
rtRvan fc Norris Ill99lh and -1.70 dex ([Carney et al. I 
119961 1. However, it is more metal poor than the halo 
metallicity derived using the SDSS F and G stars (Al- 
lende Prieto et al. 2006; Ivezic et al. 2008). Allende 
Prieto et al. (2006) found that the halo stars exhibit a 
broad range of iro n abundances with a peak at [Fe/H] 
= —1.40, whereas llvezic et al. I l|2008h find a peak at 
— 1.46 dex with a standard deviation of 0.30 dex. Re- 
stricting our sample to z < 7 kpc, as in the Ivezic et al. 
study, we find that the mean metallicity for the halo 
stars within z — 7 kpc is [Fe/H] = —1.42. Therefore, 
our metallicity estimates for both the thick disk and 
halo are consistent with the results from the SDSS 
analysis of stars within z = 7 kpc. 

Since our results are consistent with the results 
from the SDSS, and the CFHTLS data goes deeper 
than the SDSS, we can use the CFHTLS data to 
study the stars more than 7 kpc away from the Galac- 
tic plane. Our dataset is complete up to z = 17 kpc for 



1 This sample may be contaminated by metal-rich stars 
that belong to the high velocity tail of the thin disk distri- 
bution. 



g — r — 0.4 stars and up to 11 kpc for the entire sample 
(see Section 2) . Figure 6 shows that the mean metallic- 
ity for the halo decreases with increasing distance from 
the plane of the Galaxy. However, the metallicities are 
also consistent within the errors. If the decline in halo 
metallicity is real, this trend in metallicity and the ob- 
served broad distribution may be due to the inclusion 
of an increasing number of very metal poor stars in the 
halo. 

Using spectroscopic and kinematic data from the 
SDSS Data Release 5, ICarollo et ai~l (120071 ) found a 
dichotomy in the metallicity distribution of the halo 
stars with two peaks at [Fe/H] w —1.6 and [Fe/H] 
< —2.0. They found that stars that reach large dis- 
tances in their orbits above the Galactic plane show 
a metallicity peak around —2.1 dex, and they identi- 
fied those stars as members of the outer halo. Since 
the number of these outer halo stars with retrograde 
orbits is small compared to the inner halo stars, it 
is impossible to identify them without kinematic in- 
formation. However, a comparison of the observed 
metallicity distribution for our sample of stars and 
the fitted double gaussians in Figure 5 shows that the 
metal— poor tail end of the gaussian distribution for 
the halo underestimates the number of observed stars 
for z > 4 kpc. The majority of the panels in Figure 5 
show a nearly flat distribution of stars between [Fc/H] 
= —1.70 and —2.0. This could be caused by a con- 
tribution from the outer halo stars. However, without 
kinematic information, we cannot confirm or reject this 
idea. In addition, the systematic effects for the pho- 
tometric metallicit y estimates are large st for the most 
metal— poor stars, llvezic et al. I IMPOST ) find that the 
root-mean-square scatter of the residuals between the 
photometric and spectroscopic metallicity estimates of 
the F and G stars in the SDSS increases up to 0.30 dex 
for low metallicity stars. Even though we corrected for 
the differences between the MegaCam u filter and the 
standard filter set, a small error in this transformation 
may also result in systematic offsets. For example, an 
offset of 0.04 mag in the u— band data can cause a 
metallicity bias of up to 0.20 dex (Ivezic et al. 2008). 
Hence, the observed distribution may be affected by 
these systematic problems. 



5 Conclusions 

We have performed a morphological selection of point 
sources and analyzed the metallicity distribution of 
distant F and G type stars in the CFHTLS D4 field. 
We have used MegaCam photometry transformed to 
the standard ugriz filter system, and used the pho- 
tometric parallax and metallicity relations from the 
SDSS to obtain metallicities and distances for our 
sample of 1315 F and G type stars. 

Li miting our study to z = 7 kpc, the distance limit 
used bv llvezic et al. I (|2008l ). we find that our estimates 
for the mean metallicity of the thick disk ([Fe/H] = 
—0.77) and halo (—1.42 dex) are consistent with the 
results from the SDSS Data Release 6 photometry 
(-0.80 dex and -1.46 dex). Since the CFHTLS data 
is deeper than the SDSS, we are also able to check 
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Table 1: Metal abundances for the thick disk and halo stars as a function of distance from the Galactic 
plane. 
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0.68 


(20.0, 27.0] 


23.5 






-1.65 


0.62 



for trends in the thick disk and halo metallicity distri- 
bution beyond 7 kpc. Unlike the Ivezic et al. study, 
we do not find a trend in metallicity for the thick disk 
between 1 and 4 kpc. However, the significance of 
the trend observed in the SDSS dataset comes from 
the stars within 1 kpc of the Galactic plane. Unfor- 
tunately, those stars are saturated in the CFHTLS 
observations, and the metallicity trend for the thick 
disk stars further away than 1 kpc from the plane is 
consistent with the SDSS results. 

We observe a decline in metallicity of the halo with 
increasing distance from the Galactic plane. We sug- 
gest that this decline may be due to an increasing con- 
tribution from the very metal poor ([Fe/H]< —2) halo 
stars, but we cannot rule out systematic problems in 
our analysis. Proper motion data from the CFHTLS 
observations will be available at the end of the survey 
period. These data will be useful to study the kine- 
matic properties of our sample of stars and constrain 
the contribution of very metal poor stars. 
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